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\\  c  have  investigated  the  effects  of  spatial  and  temporal  variations.in  the  he  mng  of  cool-type  solar  loop 
models.  In  contrast  to  hot  loops,  these  cool  loops  have  temperature  gradients  that  are  everywhere  shallow  and 
maximum  temperatures  that  are  well  below  10"  K  The  results  of  our  numerical  simulations  show  that  spatial 
asymmetries  in  the  energy  input  to  cool  loops  will  produce  steady  state,  end-to-end  flows  Velocities  exceeding 
20  km  s  *  are  possible.  In  no  ease,  however,  is  the  downflowing  leg  appreciably  faster  or  appreciably  brighter 
than  the  upflowing  leg  tut  a  given  temperature),  and  no  significant  Dopper  shift  would  be  seen  from  an  unre¬ 
solved  loop.  Other  simulations  show  that  symmetric,  but  sudden,  reductions  in  the  energy  input  can  also 
produce  downflows  of  large  velocity.  These  flows  are  transient,  and  because  thev  arise  from  cooling  of  the 
plasma,  the  greatest  velocities  t  >4  km  s  ')  are  restricted  to  temperatures  below  5  x  I04  K.  We  conclude  from 
these  results  that  spatial  and  temporal  heating  variations  of  the  types  considered  here  are  unable  to  cvlam 
the  net  redshifts  observed  in  emission  lines  formed  near  !0'  K  on  the  Sun. 
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I.  IN  tRODl  t  rios 

\  variety  of  rocket  and  satellite  missions  have  now  estab¬ 
lished  the  existence  of  a  gener  I  downflow  into  both  quiet  and 
active  areas  on  the  Sun  Ultraviolet  emission  lines  formed  at 
temperatures  near  10'  k  in  the  transition  region  exhibit  disk- 
averaged  redshifts  of  '7  km  s  1  in  the  quiet  Sun  (Dere. 
Bartoe.  and  Brueekner  19X4.  19X6;  Roussel-Dupre  and  Shine 
19X2.  Doschek.  Feldman,  and  Bohlin  1976;  Shine  19X51  and 
-  12  km  s  1  in  active  regions  (Feldman.  Cohen,  and  Doschek 
19X2;  Dere  19X2;  Brueekner  19X1;  Klimchuk  19X6.  19X7).  Asa 
rule,  the  redshifts  are  greater  where  the  magnetic  held  is 
enhanced  (eg.,  the  network)  and  in  observations  made  closer 
to  disk  center.  We  thus  expect  the  vertical  component  of  dow  n¬ 
flow  m  magnetic  regions  to  be  somewhat  greater  than  the  disk- 
averaged  redshifts. 

Several  different  suggestions  have  been  made  as  to  the  origin 
of  (he  observed  redshifts.  One  idea  involves  the  injection  of 
cool  (  -  1()4  ki  spieular  material  into  the  corona,  which  then 
heats  and  falls  to  the  surface  at  a  higher  temperature  ( Pneuman 
and  Kopp  I97X;  Athay  and  Hol/er  19X2).  Indeed,  crude  esti¬ 
mates  of  the  global  mass  fluxes  in  the  spieular  and  UV- 
emimng  material  seem  to  agree.  No  strong  conclusions  can  be 
drawn  from  this  apparent  agreement,  however,  since  the  flux 
estimates  are  each  uncertain  by  a  factor  of  at  least  3. 

f  )t her  ideas  have  considered  the  ellects  of  spatial  and  tempo¬ 
ral  variations  in  the  healing  of  existing  hot  material  within 
magnetically  closed  regions  of  the  atmosphere  Boris  and 
Mariska  (19X2).  Mariska  ami  Boris  (19X3).  Craig  and  MeCly- 
montll9X6).  McClymont  and  Craig  ( 19X6.  19X7).  and  Mariska 
1 19X76.  {  |  studied  coronal  loop  models  and  found  that  steady 
state,  end-to-end  (lows  will  occur  if  the  energy  input  is  distrib¬ 
uted  asymmetrically  about  the  loop  apex.  While  (hese  models 
can  reproduce  many  aspects  of  the  observations,  including  the 
net  redshifts.  they  t end  lo  require  loop  lengths  that  tire  verv 
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short  (  -  10'  km)  or  pressures  that  are  smaller  Ilian  observed 
( s-  10  :  dyn  cm 

Mariska  ct  at.  (19X2)  and  Mariska  |l9X7<t|  found  that  sym¬ 
metric.  but  sudden,  reductions  in  the  energy  input  to  hot  loops 
can  also  produce  downflows  of  the  correct  magnitude.  The 
problem  here  is  that  the  emission  at  10'  K  is  much  too  faint  to 
be  important.  Thus,  both  scenarios  involving  heating  varia¬ 
tions  in  hot  loops  fall  short  in  explaining  all  aspects  of  the  solar 
observations. 

Another  class  of  loop  model,  called  cool  loops,  has  recently 
been  emphasized  in  the  literature  (Antioehos  and  Noci  19X6; 
Woods.  Holzer.  and  MacGregor  19X7u:  McClymont  19X7; 
Martens  and  Kuin  19X2;  Hood  and  Priest  1979).  Temperature 
gradients  are  everywhere  shallow  in  cool  loops,  and.  depending 
on  the  functional  forms  of  the  energy  input  and  radiative 
losses,  the  maximum  temperature  is  of  order  It)'  K.  By  com¬ 
parison.  hot  loops  have  steep  temperature  gradients  in  the 
transition  region  and  reach  a  maximum  temperature  of  order 
10"  K. 

Because  of  their  temperature  structure,  cool  loops  emit 
much  more  strongly  at  temperatures  near  10'  K  than  do  hot 
loops.  It  is  speculated  that  they  may  be  an  important  com¬ 
ponent  of  the  solar  atmosphere  and  that  they  may  explain  the 
lower  temperature,  negative  slope  end  of  the  observed  emission 
measure  curve,  w  here  hot  loop  models  are  substantially  inade¬ 
quate  (Antioehos  and  Noci  19X6;  Dowdy.  Rabin,  and  Moore 
19X6).  Here  we  explore  the  hypothesis  that  cool  loops  can  also 
explain  the  observed  redshifts. 

Because  cool  loops  are  a  relatively  new  topic  in  solar 
physics,  their  dynamics  is  only  now  being  seriously  investi¬ 
gated  (Klimchuk  and  Mariska  19X7;  McClymont  19X7).  As  an 
obvious  extension  of  the  hoi  loop  modeling  described  above, 
we  have  examined  the  dependence  of  cool  loop  models  on  both 
sptii ial  and  temporal  variations  in  the  energy  input.  Specifi¬ 
cally,  we  have  simulated  the  response  of  initially  static  cool 
loops  (o  changes  in  the  energy  input.  As  for  hot  loops,  it  is 
found  that  spatially  asymmetric  changes  produce  a  final  steady 
state  that  is  dynamic  and  that  spatially  symmetric  changes 
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produce  a  final  state  that  is  static.  We  discuss  these  two  cases 
separately  in  $$  IV  and  V.  The  numerical  model  used  for  this 
investigation  is  described  in  ij  II.  In  ij  III  we  discuss  some 
general  properties  of  cool  loop  equilibria,  paying  particular 
attention  to  the  relationship  between  structure  and  energy 
input.  Finally,  in  5;  VI  we  review  and  comment  on  our  results. 

it.  rut-  MODI  ! 

The  numerical  model  we  use  to  study  the  behavior  of  cool 
loop  atmospheres  was  developed  at  the  Naval  Research 
Laboratory  initially  for  studying  hot  loops  (e.g.,  Mariska  el  ul. 
1982;  Klimchuk.  Antiochos.  and  Mariska  1987).  The  computer 
code  solves  the  one-dimensional,  time-dependent  equations  of 
mass,  momentum,  and  energy  transport  for  a  fully  ionized 
solar  plasma  confined  within  a  rigid  magnetic  flux  tube: 

c/i  <' 

=  0  .  (1) 

<7  IS 

c  ('  ,  ('  P 
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Here.  I  is  the  time;  s  is  the  spatial  coordinate  along  the  axis 
of  the  tube;  />.  P.  T,  and  t  are  the  total  mass  density,  total 
pressure,  temperature,  and  bulk  velocity,  respectively;  ne  and 

are  the  number  densities  of  electrons  and  protons;  gv  is  the 
gravitational  acceleration  parallel  to  the  tube;  k0  is  the  coeffi¬ 
cient  of  thermal  conductivity;  and  y  is  the  ratio  of  specific 
heats,  taken  to  be  5  3  throughout  the  atmosphere.  The  term 
involving  A(T)  describes  the  rate  of  energy  loss  due  to  radi¬ 
ation.  and  e  describes  the  rate  of  energy  input  due  to  mechani¬ 
cal  heating  of  unknown  origin.  For  these  calculations  we 
assume  that  the  energy  input  depends  only  on  position, 
density,  and  time.  We  further  assume  that  the  heating  mecha¬ 
nism  does  not  impart  momentum  directly  to  the  plasma, 
though  this  may  be  a  poor  approximation  for  certain  forms  of 
wave  heating  (see  §  VI). 

The  system  of  equations  is  completed  with  the  ideal  gas  law, 

P  =  nkT  .  (5) 

where  u  is  the  total  number  density  and  k  is  Boltzmann’s  con¬ 
stant. 

We  solve  these  equations  on  a  finite-difference  grid  using 
time-step  splitting.  Convective  terms  are  treated  explicitly 
using  a  flux-corrected  transport  (FCT)  algorithm  (Book,  Boris, 
and  Hain  1975),  and  thermal  conduction  is  treated  implicitly  to 
maintain  numerical  stability.  Because  of  the  shallow  gradients 
characteristic  of  cool  loop  atmospheres,  a  computational  cell 
size  ranging  from  100  km  in  the  chromosphere  to  10  km  in  the 
overlying  regions  provides  adequate  spatial  resolution. 

The  flux  tube  geometry  is  that  of  a  semicircular  loop 
attached  to  long,  straight  chromospheric  legs.  For  the  short 
loop  we  consider  first  the  loop  radius,  and  therefore  its  height. 
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is  1.5  x  lfr'  km.  A  taller  loop  of  radius  6  x  10-'  km  is  con¬ 
sidered  in  §  Vh.  The  chromospheric  legs  range  in  size  from 
4  x  ID3  km  for  the  uniform  heating  simulations  to  8  x  I03  km 
for  the  asymmetric  heating  simulations.  In  all  cases  the  cross- 
sectional  area  of  the  flux  tube  is  constant. 

The  detailed  radiative  loss  properties  of  the  lower  transition 
region  are  not  well  understood  at  this  time.  While  it  is  clear 
that  the  atmosphere  below  '6  x  104  K  is  not  optically  thin  in 
some  transitions,  the  degree  to  which  optical  depth  effects  are 
important  is  not  obvious.  We  choose  to  follow  the  suggestion 
of  McClymont  and  Canfield  (1983)  that  the  radiative  loss  curve 
between  2  x  104  and  105  K  obeys  the  approximate  power  law 

A(T)  =  A„T3.  (6) 

A  value  for  A0  of  6.46  x  10  37  ergs  cm3  s"  1  K  '  3  gives  agree¬ 
ment  with  Raymond's  (1979)  optically  thin  loss  curve  at  103  K. 
Athay  (1986u)  has  recently  argued  that  the  optical  depth  effects 
implicit  in  equation  (6)  may  be  overemphasized  and  that  the 
T3  dependence  may  be  too  strong.  In  any  event,  it  is  possible,  if 
not  probable,  that  the  dependence  is  stronger  than  T2  above 
3  x  I04K. 

The  detailed  temperature  structure  of  the  solar  chromo¬ 
sphere  is  determined  almost  entirely  bv  the  radiative  properties 
of  the  chromospheric  plasma  (Athay  1985,  19866).  A  plateau 
occurs  near  6500  K  due  to  the  strong  temperature  dependence 
of  hydrogen  ionization  and  the  availability  of  free  electrons  to 
excite  Ca  it  and  Mg  n  transitions  responsible  for  most  of  the 
radiation.  As  a  result,  the  chromosphere  behaves  something 
like  a  thermostat,  maintaining  uniform  temperatures  over  a 
wide  range  of  pressures  and  energy  inputs.  A  second,  narrower 
temperature  plateau  occurs  near  2  x  104  K  due  to  the  escape 
of  Lyz  photons. 

For  our  model  we  have  adopted  an  isothermal  chromo¬ 
sphere  corresponding  to  a  single  temperature  plateau.  As  in 
our  earlier  hot  loop  simulations  (Klimchuk.  Antiochos,  and 
Mariska  1987)  we  assume  a  value  of  2  x  104  K.  corresponding 
to  the  upper  plateau.  In  retrospect,  we  see  that  a  choice  of  6500 
K  may  have  been  more  appropriate,  since  the  dynamical 
response  of  the  chromosphere  is  determined  mostly  by  the 
gravitational  scale  height  of  the  lower,  much  thicker  plateau. 
Our  basic  conclusions  should  not  be  adversely  affected, 
however  (2  x  104  K  remains  the  best  choice  for  hot  loop 
models,  since  for  these  models  thermal  conduction  is  impor¬ 
tant  and  should  approximately  vanish  at  the  upper  plateau). 
As  discussed  in  Klimchuk,  Antiochos,  and  Mariska  (1987).  the 
model  chromosphere  is  kept  nearly  isothermal  by  introducing 
a  precipitous  decrease  in  the  radiative  loss  curve  for  tem¬ 
peratures  below  2  x  104  K;  whenever  the  plasma  cools  to  a 
lower  temperature,  the  radiative  losses  are  severely  suppressed 
and  the  atmosphere  reheats. 

III.  STATIC'  FQUII.IBRIA 


a)  Basic  Properties 

It  is  worthwhile  to  now  review  some  of  the  basic  properties 
of  static  cool  loop  equilibria.  Further  details  can  be  found  in 
Antiochos  and  Noci  (1986)  and  in  Woods,  Holzer,  and  Mac¬ 
Gregor  (1987a),  We  first  consider  a  loop  that  is  heated  sym¬ 
metrically  by  an  energy  input  that  is  uniform  in  space.  Figure  1 
shows  the  equilibrium  profiles  of  temperature,  density,  and 
pressure  that  result  for  a  quict-Sun  base  pressure  (at  the  top  of 
the  chromosphere)  of  0.2  dyn  cm  2  (Mariska  1986).  The  apex 
temperature  is  7.5 
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Fui.  i.  -  The  static  equilibrium  structure  of  a  uniformly  heated  short  loop  of  radius  1.5  x  !03  km:  temperature  (K),  total  pressure  (dyn  cm  ~2),  and  total  number 
density  (cm  t  are  plotted  as  a  function  of  position  along  the  loop  axis.  Only  half  of  the  loop  is  shown  because  of  symmetry.  The  apex  is  located  at  position 
1.04  x  104  km  and  the  base  (top  of  chromosphere)  at  position  8.0  x  103  km.  Please  note  that  the  vertical  scale  is  logarithmic. 


Unlike  hot  loops,  whose  temperatures  rise  first  abruptly  in 
the  transition  region  and  then  very  minimally  in  the  corona, 
static  cool  loops  have  temperatures  that  increase  both  grad¬ 
ually  and  steadily  throughout.  This  contrast  is  related  to  the 
fundamental  difference  in  energy  balance  in  these  two  types  of 
atmospheres.  In  hot  loops,  thermal  conduction  plays  an 
important  role  transporting  energy  from  the  corona  down  to 
the  transition  region;  steep  temperature  gradients  are  neces¬ 
sary  at  the  lower  temperatures  to  maintain  the  large  conduc¬ 
tive  flux.  In  cool  loops,  thermal  conduction  is  relatively 
unimportant,  and  there  is  everywhere  a  local  balance  between 
radiation  and  energy  input. 

In  a  uniformly  heated  cool  loop,  radiation  is  also  uniform, 
and  the  temperature  increases  with  height  in  order  to  offset  the 
effects  of  decreasing,  gravitationally  stratified  density 
[n2A(T)  =  constant].  If  the  heating  decreases  with  height,  the 
temperature  rises  less  quickly,  and  the  heating  increases  with 
height,  the  temperature  rises  more  quickly.  For  radiative  loss 
curves  having  a  weaker  dependence  than  T2,  the  heating  must 
not  only  decrease  with  height,  but  it  must  decrease  sufficiently 
rapidly;  otherwise,  cool  equilibria  are  not  possible  (Woods, 
Holzcr,  and  MacGregor  1987a). 

The  maximum  temperature  attainable  by  a  cool  equilibrium 
under  the  conditions  of  our  model  is  105  K,  the  temperature  at 
which  the  constant  pressure  form  of  the  loss  curve,  A(T)/T2, 
first  has  a  negative  slope.  Beyond  this  a  temperature  increase 
cannot  offset  a  density  decrease,  and  the  radiation  is  necessar¬ 
ily  reduced.  Energy  balance  is  impossible  with  uniform  heating 
and  it  is  difficult  with  heating  that  decreases  with  height,  since 
AIT)  is  itself  negatively  sloped  beyond  10s  K. 

The  maximum  temperature  in  our  model  places  an  upper 
limit  on  the  possible  height  of  the  loop.  With  heating  that  is 
uniform  in  space,  the  maximum  height  is  2  x  103  km,  and  with 
heating  that  is  constant  per  unit  mass,  the  maximum  height  is 
8  x  I03  km  (Antiochos  and  Noci  1986).  Greater  temperatures 
and  heights  are  possible  for  different  loss  curves  and  heating 


functions,  however  (Woods,  Holzer,  and  MacGregor  1978 a). 
There  is  no  limit  to  the  height  of  hot  loops. 

The  gravitational  scale  height  of  a  hot  loop  is  generally  large 
compared  to  its  geometric  height,  and  consequently  the  pres¬ 
sure  is  nearly  constant.  The  same  is  not  true  of  cool  loops,  on 
the  other  hand,  which  have  a  noticeable  pressure  stratification. 
In  Figure  1,  for  example,  the  pressure  decreases  by  a  factor  of  2 
from  the  base  to  the  apex.  The  density  decreases  even  more 
because  of  the  temperature  increase  required  for  energy 
balance. 

b)  Dependence  on  Energy  Input 
We  now  examine  the  relationship  between  atmospheric 
structure  and  energy  input  in  more  detail.  It  will  be  assumed 
here  that  the  energy  input  depends  only  on  position. 

Thermal  conduction  is  fully  ignorable  for  the  cool  loops  we 
are  now  considering,  and  the  static  energy  equation  reduces  to 

A0n2T3  =  «  ,  (7) 

where  we  have  taken  ne  =  np  =  n/2.  With  the  aid  of  the  ideal 
gas  law,  this  can  be  rewritten  as 


showing  that  temperature  is  directly  proportional  to  energy 
input.  In  regions  of  uniform  pressure,  the  temperature  and 
energy  input  profiles  will  therefore  have  similar  form.  Local 
deviations  in  heating  should  produce  local  deviations  in  tem¬ 
perature. 

The  static  momentum  equation  is  simply 
dP 

—  =  -  mpneg  ,  (9) 

where  mp  is  the  proton  rest  mass,  z  is  the  vertical  spatial  coor¬ 
dinate.  increasing  upward,  and  g  is  a  positive  constant.  The  z 
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coordinate  should  not  be  confused  with  the  curvilinear  coordi¬ 
nate.  s,  in  equations  { 1 )— (3).  Substituting  for  nr  using  the  ideal 
gas  law  and  equation  (8)  we  get 


dP  mpcjA0  P3 

7^=  8  k3e 


HO) 


The  pressure  at  any  height,  including  the  apex,  h,  can  then  be 
easily  obtained  by  integrating  upward  from  the  base  (top  of 
chromosphere): 


P[h)  =  P  o 


-  1  2 

(11) 


where  PQ  is  the  base  pressure. 


(12) 


amplitude  equal  to  50%  of  the  background  heating  rate.  The 
source  is  activated  over  a  period  of  100  s,  after  which  time  the 
asymmetric  heating  is  held  constant.  Figure  2  indicates  the 
response  of  the  loop  to  the  additional  heating:  temperature 
and  velocity  profiles  are  shown  for  selected  times  both  before 
and  after  the  source  is  turned  on.  (Note  that  this  and  following 
temperature  plots  in  this  paper  are  plotted  linearly,  rather  than 
logarithmically,  as  is  usually  done  for  hot  loops.) 

Temperatures  rise  quickly  in  the  region  of  the  additional 
heating,  and  within  300  s  a  bump  appears  on  the  temperature 
profile  having  a  similar  shape  and  magnitude  to  the  energy 
source  itself.  This  is  expected  from  equation  (8),  since  pressure 
is  roughly  uniform  over  the  region  and  the  flows  are  not  large 
enough  to  significantly  affect  the  equilibrium  structure.  Else¬ 
where  in  the  loop  the  temperatures  are  essentially  unchanged, 
reflecting  the  fact  that  the  energy  input  in  those  regions  is  also 
unchanged. 

Mass  motions  produced  by  the  additional  heating  are  ini- 


is  the  gravitational  scale  height  at  the  base,  and 

cl:) 

/<-')  =  —  (13) 

«o 

is  the  heating  function  normalized  to  its  base  value.  Equation 
( 1 1 )  is  a  simplified  form  of  the  more  general  equation  derived 
by  Antiochos  and  Noci  (1986,  their  eq.  [17]). 

We  see  from  equation  (11)  that  the  pressure  at  any  given 
height  in  a  cool  loop  depends  on  the  heating  at  all  levels  below, 
but  not  on  the  heating  above.  It  also  depends  on  the  base 
pressure,  which  is  given  uniquely  by  the  base  heating  (and  base 
temperature)  according  to  equation  (8).  This  result  contrasts 
with  hot  loops,  where,  because  of  thermal  conduction,  the  pres¬ 
sure  at  each  height  depends  on  the  heating  at  all  other  loca¬ 
tions  throughout  the  loop. 

In  order  for  a  cool  loop  to  be  in  static  equilibrium,  the  apex 
pressure  given  by  equation  (11)  must  be  the  same  for  both  of 
the  loop  legs.  Otherwise,  an  unbalanced  pressure  gradient 
would  exist  at  the  top,  where  gravity  exerts  no  force  parallel  to 
the  loop  axis.  Assuming  equal  base  pressure  at  the  two  foot- 
points,  pressure  equilibrium  at  the  apex  is  trivially  satisfied  in 
the  case  of  symmetric  heating.  For  asymmetric  heating,  on  the 
other  hand,  a  pressure  imbalance  will  drive  a  flow  from  the 
higher  pressure,  more  strongly  heated,  leg  to  the  lower  pres¬ 
sure,  less  strongly  heated,  leg.  Any  equilibrium  must  necessar¬ 
ily  be  dynamic. 

IV.  ASYMMETRIC  HEATING 

In  principle,  one  could  study  the  dynamic  equilibria  that 
result  from  asymmetric  heating  by  solving  the  steady  state 
equations  directly.  Instead,  we  obtain  the  steady  state  solution 
by  applying  asymmetric  heating  to  an  initially  static  cool  loop 
and  allowing  it  to  evolve  asymptotically  to  a  final  equilibrium. 
In  that  way  we  can  also  examine  the  intermediate,  nonequilib¬ 
rium  states  that  might  be  expected  to  occur  on  the  Sun. 

a /  Localized  F.nergv  Source 

In  the  first  of  our  simulations  we  consider  the  efieets  ,'f 
adding  energy  locally  to  only  one  side  of  a  loop.  Beginning 
with  the  static  atmosphere  of  Figure  1,  and  without  altering  the 
uniform  background  heating,  we  gradually  turn  on  an  energy 
source  concentrated  haifway  up  one  of  the  iegj.  ».  iia.»  a  Gauss¬ 
ian  distribution  with  an  e  1  half-width  of  500  km  and  a  central 


tially  away  from  the  source  in  both  directions.  This  is  only  a 
short  term  response  to  the  local  pressure  increase,  however, 
and  within  200  s  the  atmosphere  adjusts  to  set  up  a  unidirec¬ 
tional  flow  away  from  the  footpoint  on  the  more  strongly 
heated  side.  Flow  velocities  gradually  increase  until  an  approx¬ 
imate  steady  state  is  achieved  at  t  =  1000  s.  After  this  time, 
until  the  simulation  is  terminated  at  f  =  7000  s,  the  atmosphere 
is  virtually  constant.  A  maximum  velocity  of  5  km  s~  1  occurs 
in  the  upflowing  leg.  at  the  approximate  location  of  the  energy 
source.  The  striking  similarity  between  the  equilibrium  velocity 
and  temperature  profiles  is  a  direct  consequence  of  mass  con¬ 
servation  in  a  region  of  slowly  varying  pressure. 

The  flow  produced  here  is  too  slow  to  explain  the  net  red- 
shifts  observed  on  the  Sun,  regardless  of  the  emission  proper¬ 
ties  of  the  plasma.  McCIymont  and  Craig  (1987)  and  Mariska 
(19876,  c )  have  recently  shown  that  in  asymmetrically  heated 
hot  loops  the  flow  velocities  are  greatest  when  the  energy  input 
is  highly  concentrated  in  a  narrow  region  close  to  the  chromo¬ 
sphere.  To  determine  whether  this  is  also  true  of  cool  loops  we 
have  performed  a  second  simulation,  similar  to  the  first,  but 
with  the  Gaussian  energy  source  centered  only  300  km  above 
the  base  and  having  a  half-width  of  only  100  km.  The  central 
amplitude  is  200%  of  the  background  heating  rate,  which  is 
close  to  the  upper  limit  allowed  by  the  cool  state  (for  greater 
amplitudes  the  temperature  rises  above  105  K  and  the  atmo¬ 
sphere  quickly  evolves  to  a  I06  K.  hot  state). 

The  effects  of  this  type  of  heating  are  qualitatively  similar  to 
before.  Temperatures  become  enhanced  at  the  source,  this  time 
forming  a  spike  on  the  profile  rather  than  a  bump,  and  a  steady 
state  flow  is  produced  carrying  material  away  from  the  more 
strongly  heated  leg.  At  no  time  do  the  velocities  exceed  3  km 
s'1,  however.  We  conclude  that  neither  the  location  nor  the 
concentration  of  the  energy  source  is  a  strong  factor  in  deter¬ 
mining  the  maximum  flow  speed  in  cool  loops. 

b)  Exponential  Form 

In  both  of  these  simulations  the  energy  content  of  the  source 
is  small  compared  to  the  spatially  integrated  background 
heating  of  the  loop.  Consequently,  the  left  leg  receives  only 
~  20%  more  total  energy  than  the  right  leg.  To  determine 
whether  an  asymmetry  in  the  total  amount  of  heating  is  impor¬ 
tant  for  driving  flows,  we  have  considered  a  completely  differ¬ 
ent  form  for  the  heating  profile.  Instead  of  a  localized  energy 
source,  we  introduce  a  component  ot  healing  that  uc>,rcaaCj 
exponentially  along  the  axis  of  the  loop,  from  one  footpoint  to 
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Position  (103  km) 

Fio  2.  The  response  of  the  static  atmosphere  of  Fig.  I  to  an  imposed  heating  asymmetry :  profiles  of  temperature  (a)  and  velocity  (b)  are  shown  for  selected  times 
after  a  localized  energy  source  has  been  turned  on  at  position  9.2  x  103  km  (half-way  up  the  left  leg).  Positive  velocities  are  in  the  direction  of  increasing  position 
lupflow  in  the  left  leg  and  downflow  in  the  right).  The  temperature  profiles  correspond  to  times  ofO,  100.  and  1000  s.  and  the  velocity  profiles  correspond  to  times  of  0. 


100.  200.  500.  and  I (XX)  s.  After  1000  s  the  atmosphere  is  essentially  unchanged.  i.< 


the  other.  For  a  characteristic  damping  length  equal  to  one- 
half  of  the  loop  length  and  a  magnitude  equal  to  1.5  times  the 
background  heating  rate  at  the  left  footpoint,  the  left  leg 
receives  nearly  twice  as  much  total  energy  as  the  right. 
(Substantially  larger  differences  are  not  possible  without  the 
atmosphere  evolving  to  a  hot  state.) 

As  we  anticipate  on  the  basis  of  equation  (II),  the  flow  pro¬ 
duced  by  this  heating  asymmetry  is  considerably  faster  than  in 
the  previous  two  cases.  Figure  3  shows  the  beginning,  interme¬ 
diate.  and  steady  state  profiles  of  temperature,  velocity,  and 
total  number  density.  The  profiles  are  less  structured  than  in 
the  previous  simulation';  because  the  energy  input  varies  inure 
slowly  with  position.  Also,  the  maxima  in  temperature  and 
velocity  now  occur  on  the  side  of  the  loop  with  lesser  heating. 


steady  state. 


rather  than  on  the  side  with  greater  heating.  This  occurs 
because  the  heating  increases  with  height  in  the  right  leg,  but 
decreases  with  height  in  the  left ;  temperature  gradients  must  be 
steeper  in  the  right  leg  so  that  radiation  can  balance  the  energy 
input  (§  Ilia). 

The  chromosphere  of  the  final  equilibrium  is  depressed  at 
both  the  left  and  right  footpoints  compared  to  the  initial  state. 
This  is  understood  in  terms  of  the  energy  balance  equation, 
equation  (8),  with  e  taken  to  be  the  effective  heating  rate, 
including  the  contribution  due  to  flow  terms.  Comparing  the 
initial  and  final  states  we  find  that  the  effective  heating  rate  is 
increased  at  both  of  the  footpoims.  At  the  right  footpoint,  flow 
heating  augments  the  small  increase  in  mechanical  heating, 
and  at  the  left  footpoint,  flow  cooling  is  insufficient  to  offset  the 
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Fig  3.— The  response  of  the  static  atmosphere  of  Fig.  1  to  an  imposed 
heating  asymmetry  profiles  of  temperature  (a),  velocity  (b),  and  total  number 
density  |c)  are  shown  for  selected  limes  after  an  energy  source  has  been  turned 
on  that  decreases  exponentially  from  the  left  footpoint  to  the  right  along  the 
axis  of  the  loop.  The  temperature  and  density  profiles  correspond  to  times  of  0, 
100  {dotted  line),  and  1 000s,  and  the  velocity  profiles  correspond  to  times  of  0, 
100.  300.  500,  and  1000  v  After  1000  s  the  atmosphere  is  essentially  unchanged, 
i  e..  steady  state  Note  that  the  density  scale  is  logarithmic. 

large  increase  in  mechanical  heating.  On  both  sides  the  magni¬ 
tude  of  the  flow  term  is  roughly  50%  of  the  original  mechanical 
heating. 

Noting  that  changes  in  base  height  are  related  to  changes  in 


base  pressure  according  to  A z0/H0  -AP0/P0,  we  see  from 

equation  (8)  that 


At0  H  o 
*0  2  • 


(14) 


where  €0  is  the  effective  heating  and  H0  is  the  gravitational 
scale  height.  For  Ac ,  positive,  Ar0  is  negative  and  the  chromo¬ 
sphere  is  depressed,  as  observed  in  the  simulation. 

The  large  velocities  of  order  20  km  s  '1  in  Figure  3  suggest 
that  the  net  redshil'ts  observed  on  the  Sun  may  be  a  result  of 
asymmetric  heating  of  this  type.  Both  downflow  and  upflow 
are  present,  however,  and  a  net  redshift  will  be  produced  only  if 
the  downflowing  leg  is  appreciably  faster  or  appreciably 
brighter  than  the  upflowing  leg.  Comparing  the  velocity  and 
temperature  profiles  of  Figure  3,  we  see  that  the  downward  and 
upward  flow  speeds  are  similar  at  each  temperature.  A  net 
redshift  is  possible,  therefore,  only  if  the  downflowing  leg  is 
relatively  bright. 

The  brightness  of  an  optically  thin  emission  line  is  pro¬ 
portional  to  the  emission  measure,  or  the  square  of  the  electron 
density  integrated  over  the  emitting  volume.  This  emitting 
volume  is,  in  turn,  related  to  the  inverse  temperature  gradient 
through  the  strong  dependence  of  ionization  on  temperature. 
As  Figure  3  indicates,  the  steady  state  densities  in  the  loop  are 
a  similar  function  of  temperature  in  both  legs.  Differences  in 
emission  measure  are  therefore  related  to  differences  in  tem¬ 
perature  gradient  alone. 

Since  temperature  gradients  are  shallower  in  the  upflowing 
leg  of  the  loop,  we  expect  that  leg  to  be  brighter.  This  is 
verified  in  Figure  4,  where  we  have  plotted  the  emission  mea¬ 
sures  calculated  for  temperature  intervals  of  0.3  in  the 
logarithm  (A  log  T  —  0.3)  at  the  time  of  maximum  velocity  in 
Figure  3;  upflowing  and  downflowing  legs  are  given  separately. 
The  upflows  are  brighter  than  the  downflows  at  all  but  the 
highest  temperatures,  where  they  are  comparable.  Thus, 
despite  the  large  velocities,  the  heating  asymmetry  considered 
here  is  unable  to  reproduce  the  net  redshifts  observed  on  the 
Sun.  If  anything,  this  simulation  would  predict  a  net  blueshift. 


V.  SYMMETRIC  HEATING  VARIATIONS 

We  now  examine  the  effects  of  heating  variations  that  are 
symmetric  about  the  loop  apex.  Whereas  asymmetric  varia¬ 
tions  produce  flows  directed  from  one  end  of  the  loop  to  the 
other,  symmetric  variations  must  produce  flows  that  are  of  the 
same  sense  in  both  legs,  either  downward  or  upward.  In 
general,  we  expect  cooling  atmospheres  to  drain,  giving  rise  to 
downflows,  and  heating  atmospheres  to  fill,  giving  rise  to 
upflows.  This  is  not  always  the  case,  however,  as  we  demon¬ 
strate  shortly. 

To  study  these  effects  we  have  performed  a  series  of  cooling 
and  heating  simulations  similar  to  those  performed  by 
Mariska  (1987a)  for  hot  loops.  Beginning  with  a  static  atmo¬ 
sphere,  we  symmetrically  reduce  the  energy  input  by  a  pre¬ 
scribed  amount  over  a  prescribed  period  of  time.  While 
maintaining  this  reduced  input  we  then  allow  the  atmosphere 
to  evolve  toward  the  new  static  equilibrium.  Once  the  flows 
have  diminished  to  well  under  1  km  s' ',  we  restore  the  input 
to  its  original  level,  and  the  atmosphere  returns  to  its  initial 
state.  In  this  way  the  loop  passes  through  a  complete  cycle  of 
cooling  and  heating.  Note  that  only  half  the  loop  need  be 
modeled  because  of  symmetry. 


Log  Temperature  (K) 

F Ki.  4  Plot  of  emission  measure  versus  temperature  near  the  time  of  maximum  velocity  1 1000  s)  in  the  asymmetric  heating  simulation  of  Fig.  3.  The  upflowing 


0  {solid  line  I  and  downflowing  (dotted  line )  legs  are  given  separately.  The  units  of  emission  measure  are  arbitrary. 


a)  Short  Loop 

The  relatively  short  loop  we  consider  first  is  identical  to  that 
used  in  the  previous  simulations:  its  radius  of  curvature,  and 
therefore  height,  is  1.5  x  103  km.  In  the  cooling  half  of  the 
cycle  the  uniform  energy  input  is  everywhere  reduced  by  90% 
over  a  period  of  100  s.  Figure  5  shows  the  atmospheric 
response 

The  atmosphere  cools  quickly,  on  a  time  scale  that  is  compa¬ 
rable  to  the  radiative  cooling  time: 


t, - -  200  s  . 

*  (y-l)A  nPT 


This  will  always  be  the  case  as  long  as  the  heating  is  reduced 
over  a  time  scale  that  is  short  compared  to  the  radiative 
cooling  time.  If  the  heating  is  instead  reduced  over  a  longer 
time  scale,  the  atmosphere  will  cool  at  the  rate  of  the  reduction 
itself. 

Initially  the  flow  is  downward  in  the  loop  as  the  atmosphere 
collapses  from  the  decrease  in  thermal  support.  The  chromo¬ 
sphere  does  not  cool,  however  (because  of  the  isothermal 
assumption  built  into  the  radiative  loss  curve),  and  it  expands 
upward  into  the  region  of  decreasing  pressure.  During  this 
expansion  the  chromosphere  acts  as  a  snowplow,  displacing 
the  hotter  material  in  its  way.  Before  long  the  entire  loop  is 
experiencing  an  upflow.  The  greatest  velocities  occur  at  the 


very  top  of  the  chromosphere,  but  at  all  times  they  remain 
slower  than  4  km  s  '. 

Within  2000  s  the  atmosphere  has  nearly  settled  into  a  new 
static  equilibrium.  We  then  restore  the  energy  input  to  its  orig¬ 
inal  level.  The  resulting  behavior  is  analogous,  but  opposite,  to 
before:  temperatures  and  pressures  rise  in  the  loop,  and  the 
heated  plasma  expands,  forcing  the  chromosphere  downward. 
After  this  initial  adjustment  a  small  upflow  develops  as  the 
loop  slowly  fills  with  material  to  regain  its  original  densities. 

The  maximum  downflow  of  the  entire  simulation  occurs 
during  the  cooling  phase  of  the  cycle.  Because  it  is  so  slow,  and 


because  it  involves  lower  temperature  material,  the  observed 
solar  redshifts  cannot  be  explained. 

h)  Tall  Loop 

For  two  reasons  we  expect  that  a  significantly  taller  loop  will 
respond  differently  to  a  similar  decrease  in  heating:  first,  the 
chromospheric  expansion,  which  is  limited  to  roughly  1  gravi¬ 
tational  scale  height  according  to  equation  (14),  will  influence 
only  the  lower  portion  of  a  tall  loop,  and  second,  the  material 
near  the  top  of  a  tall  loop  will  fall  a  greater  distance  and,  in 
principle,  should  be  accelerated  to  a  greater  velocity.  We  have 
therefore  considered  a  second  loop  that  is  4  times  larger  than 
the  first.  Its  radius  of  curvature  is  6  x  103  km. 

In  order  to  obtain  an  initial  cool  equilibrium  in  this  tall  loop 
it  is  necessary  to  assume  an  energy  input  that  decreases  with 
height;  otherwise,  only  a  hot  equilibrium  is  possible  (§  Ilia). 
Instead  of  prescribing  an  explicit  height  dependence,  we 
assume  a  heating  rate  that  is  constant  per  unit  mass  and  there¬ 
fore  decreases  implicitly  with  height: 

e{p)  =  Q0P .  ('6) 

where  Q0  is  a  constant,  set  equal  to  9.1  x  1010  ergs  s- 1  g  1  to 
achieve  a  base  pressure  of  0.2  dyn  cm  " 2. 

In  the  cooling  half  of  the  simulation  Q0  is  reduced  by  90% 
over  a  period  of  100  s.  The  loop  responds  by  draining,  as 
indicated  in  Figure  6.  Some  initial  upflow  occurs  in  front  of  the 
expanding  chromosphere,  but  it  is  quickly  overcome  by  the 
larger  scale  downflow.  A  maximum  downward  velocity  of  17 
km  s  ~ 1  is  achieved  ~  320  s  into  the  simulation. 

This  downflow  is  more  than  fast  enough  to  explain  the  solar 
redshifts;  however,  because  the  downward  acceleration  is  so 
closely  linked  to  cooling,  the  temperatures  must  drop  substan¬ 
tially  before  the  velocities  can  become  large.  This  is  made  clear 
in  Figure  7,  where  we  have  plotted  velocity  versus  temperature 
for  three  representative  times  during  the  simulation.  At  no  time 
does  material  hotter  than  5  x  104  K  exceed  a  velocity  of  5  km 
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Fig.  7.  —Plots  of  velocity  versus  temperature  at  times  of  0, 80, 160.  and  240  s  during  the  cooling  tall  loop  simulation  of  Fig,  6. 


s  Thus,  although  the  downflows  are  large,  they  are  unable 
to  explain  the  redshifts  observed  at  temperatures  near  105  K. 

An  interesting  feature  of  Figure  6  is  the  shock  that  begins  to 
form  as  the  accelerating  downflow  piles  up  material  in  the 
region  above  the  chromosphere.  The  Mach  number  is  ^0.6  at 
the  last  time  plotted  As  is  typically  the  case,  the  temperatures 
are  elevated  in  the  region  behind  the  shock  front. 

In  the  heating  phase  of  the  cycle  the  behavior  is  once  again 
analogous,  but  opposite,  to  that  of  the  cooling  phase.  The 
maximum  velocity  is  upward,  at  17  km  s' l.  Downflows  are  not 
significant. 

We  have  performed  several  other  cooling  and  heating  simu¬ 
lations  using  the  same  tail  loop,  but  varying  both  the  speed  and 
the  degree  to  which  we  vary  the  energy  input.  As  expected,  the 
flow  velocities  are  smaller  when  the  energy  input  is  varied  more 
slowly  or  when  it  is  varied  bv  a  lesser  amount.  With  the  heating 
constant,  Qn,  reduced  by  90"..  over  a  period  of  500  s.  for 
example,  the  maximum  velocity  is  only  6  km  s'1.  With  it 
reduced  by  20"..  over  100  s,  the  maximum  is  2  km  s"  *.  More 
rapid  reductions  in  the  energy  input  will  have  little  effect  on  the 
flow,  since  the  dynamics  of  the  loop  are  then  limited  by  the 
radiative  cooling  time,  which  is  of  order  200  s. 

VI.  SUMMARY  AND  DISC  USSION 

We  have  investigated  the  effects  of  spatial  and  temporal 
variations  in  the  heating  of  cool  loop  models  in  an  attempt  to 
explain  the  net  redshifts  that  are  observed  on  the  Sun.  In  none 
of  the  situations  we  have  considered  are  the  induced  flows  able 
to  satisfactorily  reproduce  the  observations.  In  the  case  of 
asymmetric  heating,  the  end-to-end  flows  can  be  as  fast  as  20 
km  s  '.  but  the  downflowing  leg  is  neither  appreciably  faster 
nor  appreciably  brighter  than  the  upflowing  leg;  no  net  red- 
shift  is  produced.  In  the  case  of  symmetric  heating,  the  down¬ 
flows  can  also  be  large,  but  they  arc  restricted  to  temperatures 
that  are  well  below  10'  K.  Neither  situation  would  give  rise  to 
the  >  7  km  s  1  disk-averaged  redshifts  seen  in  emission  lines  of 
species  like  C  tv. 

These  results  are  restricted,  of  course,  by  the  assumptions  of 


our  model.  One  major  assumption  is  the  form  of  the  radiative 
loss  curve,  which  we  have  taken  to  increase  as  T3  up  to  a 
maximum  at  105  K.  Most  of  the  optically  thin  loss  curves  that 
have  been  published  to  date  indicate  a  weaker.  -  T2.  depen¬ 
dence  between  3  x  104  and  8  x  I04  K,  although  uncertainties 
in  factors  such  as  atomic  cross  sections  and  elemental  abun¬ 
dances  give  rise  to  considerable  deviation  (Raymond  1979; 
Cox  and  Tucker  1969;  Pottasch  1965;  Summers  and  McWhir- 
ter  1979;  Gaetz  and  Salpeter  1983).  The  T3  dependence  was 
proposed  originally  by  McClymont  and  Canfield  (1983)  to 
account  for  optical  depth  effects  and  has  been  adopted  by 
other  investigators  (e.g.,  Antiochos  and  Noci  1986).  Athay 
(1986a)  has  argued  that  this  is  an  overcorrection,  and  so  the 
relevancy  of  our  results  may  be  questioned. 

To  test  the  sensitivity  of  our  results  to  the  assumed  form  of 
the  loss  curve  we  have  performed  a  limited  number  of  asym¬ 
metric  heating  simulations  using  a  loss  curve  with  a  Tl  5 
dependence.  A  background  heating  that  decreases  with  height 
is  required  for  the  initial  static  equilibrium;  we  assume  the 
density-dependent  form  of  equation  (16).  Under  these  condi¬ 
tions,  we  find  that  cool  atmospheres  are  easily  driven  to  a  hot 
state  by  even  rather  small  heating  asymmetries.  In  the  one 
cool,  steady  flow  state  that  was  achieved  there  is  not  enough 
difference  between  the  upflowing  and  downflowing  legs  to 
produce  a  significant  Doppler  shift.  Our  conclusion  that  the 
observations  cannot  be  explained  by  spatial  variations  in  the 
heating  is  therefore  independent  of  the  assumed  form  of  the 
loss  curve,  it  would  seem. 

Throughout  our  calculations  we  have  made  the  additional 
assumption  that  the  mechanism  responsible  for  the  heating 
does  not  impart  momentum  directly  to  the  atmosphere,  i.e.. 
does  not  exert  a  force  This  may  be  valid  if  the  heating  is  a 
purely  local  phenomenon  (e.g.,  conversion  of  magnetic  free 
energy);  however,  if  the  heating  is  due  to  the  dissipation  of  a 
wave  flux  then  there  is  a  force  associated  with  the  gradient  in 
the  wave  pressure.  In  their  extensive  treatment  of  the  problem. 
Woods.  Holzer.  and  MacGregor  (19876)  have  shown  that,  for 
certain  types  of  waves,  this  force  can  be  so  large  as  to  signifi- 
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cantly  affect  the  structure  of  static  equilibrium  atmospheres.  In 
fact,  such  equilibria  do  not  always  exist  when  the  wave  effects 
are  properl)  taken  into  account. 

To  see  this,  consider  the  follow  ing  heuristic  description  of  a 
wave  tluv  that  is  incident  from  below  and  damped  exponen¬ 
tially  in  the  loop: 

f(s)  =  /•,,('  'r.  (17) 

where  l.  is  the  damping  length  and  is  the  flux  at  the  base. 
The  heating  and  wan  pressure  lenergy  density  )  associated  with 
this  llux  are 


and 

P,  =  -  .  (Iff  I 

c 

respectively,  where  c  is  the  wave  speed.  The  force  exerted  by  the 
waves  is  then  given  by 


Lor  wave  speeds  that  increase  with  height,  as  in  the  case  of 
sound  waves  or  most  Alfven  waves,  the  wave  force  is  necessar¬ 
ily  upward  with  a  magnitude  of  at  least  c  c.  In  our  models 
t  '  III  :  ergs  cm  '  s  and  therefore  the  force  due  to  sound 
waves  would  be  stronger  than  the  force  due  to  gravity!  Alfven 
waves  would  likely  exert  a  weaker  force,  but  their  importance 
depends  on  the  details  of  the  particular  model. 


The  symmetric  cooling  simulations  of  t;  V7>  are  unable  to 
reproduce  the  observed  redshifts  because  in  that  model  the 
atmosphere  is  supported  solely  by  thermal  pressure  gradients: 
downflows  cannot  become  large  until  substantial  cooling  has 
occurred.  If  w  ave  pressure  effects  had  been  included,  the  initial 
atmosphere  might  have  been  supported  largely  by  these  pres¬ 
sure  gradients,  and  then  a  reduction  in  the  incident  wave  llux 
would  allow  fast  downflows  to  occur  independent  of  the 
cooling.  Thus,  there  still  remains  the  possibility  that  a  cooling 
scenario  can  explain  the  solar  observations.  We  hope  to 
explore  this  further  in  future  simulations  that  account  fully  for 
the  existence  of  waves. 
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